What is the Criterion for a Strong First Order Electroweak Phase Transition in 

Singlet Models? 
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It is widely believed that the existence of singlet scalars in some Standard Model extensions can 
easily make the electroweak phase transition strongly first order, which is needed for the electroweak 
baryogenesis scenario. In this paper, we will examine the strength of the electroweak phase transition 
in the simplest extension of the Standard Model with a real singlet using the sphaleron energy at 
the critical temperature. We find that the phase transition is stronger by adding a singlet; and also 
that the criterion for a strong phase transition Q(T C )/T C > 1, where Q — (v 2 + (x — xo) 2 )? and x 
(xo) is the singlet vev in the broken (symmetric) phase, is not valid for models containing singlets, 
even though often used in the literature. The usual condition v c /T c > 1 is more meaningful, and it 
is satisfied for a large part of the parameter space for physically allowed Higgs masses. 

PACS numbers: 98.80.Cq 11.10. Wx 11.15.Ha 



I. INTRODUCTION 

The Standard Cosmological Model has been success- 
ful in describing the early universe, it is supported by a 
number of important observations: the expansion of the 
universe, the abundance of the light elements and the 
cosmic microwave background (CMB) radiation. These 
three measurable signatures strongly support the notion 
that our universe evolved from a dense, nearly featureless 
hot gas, just as the Big Bang model predicts. 

However it fails to explain some serious problems like 
the nature of dark matter and dark energy; and the dom- 
inance of matter over antimatter. 

From a theoretical point of view, there is no justifica- 
tion to assume that the universe started its evolution with 
a defined baryon asymmetry; (t = 0) > (t = 0). 
The natural assumption is that the universe was initially 
neutral. Direct observations show that the universe con- 
tains no appreciable primordial antimatter. In addition, 
the success of big bang nucleosynthesis requires that the 
ratio of the effective baryon number (jif, — rig) to the en- 
tropy density should be between [l[ 



2.6 x 1(T 10 < j] 



rib 



< 6.2 x 10 



-10 



(1) 



This number has been independently determined to be 
rj = (8.7 ±0.3) x 10~ n from precise measurements of the 
relative heights of the first two microwave background 
(CMB) acoustic peaks by the WMAP satellite 0. Thus 
how can one understand the origin of this asymmetry? 
This is what is called the problem of baryogenesis (for 
a review see Q)- In 1967, Sakharov forwarded his three 
conditions for baryogenesis 0] , which are summarized in 
the existence of processes which: (1) violate B number, 
(2) violate C and CP symmetries; and (3) take place out 
of equilibrium. 



One of the most interesting scenarios for baryogene- 
sis is the electroweak baryogenesis (For a review, see 
Q), where the third Sakharov condition is realized via 
a strong first order phase transition at the electroweak 
scale. 

In gauge theories, a first order phase transition takes 
place if the vacuum of the theory does not correspond 
to the global minimum of the potential. Since it is en- 
ergetically unfavored, the field changes its value to the 
true vacuum (i.e. the absolute minimum of the poten- 
tial). Because of the existence of a barrier between the 
two minima, this mechanism can happen by tunneling 
or thermal fluctuations via bubble nucleation. The elec- 
troweak baryogenesis scenario is realized when the B and 
CP violating interactions pass through the bubble wall. 
These interactions are very fast outside the bubbles but 
suppressed inside. Then a net baryon asymmetry results 
inside the bubbles which are expanding and filling the 
universe at the end. 

In order to compute the net baryon number, the rate 
of B violating processes in the broken phase is needed. In 
the Standard Model, B number is violated at the quan- 
tum level where the transition between two topolog- 
ically distinct SU (2) l ground states is possible. 

The transition between two neighboring ground states 
breaks both lepton and baryon numbers by AL = AB = 
3. To find the rate of B violating processes, one needs 
to know the sphaleron solution, which is a static field 
configuration that interpolates between the two distinct 
ground states. The sphaleron configuration was found in 
H for the SU(2) L model. 

A model-independent condition in order that the phase 
transition be strong enough was derived in [f|: 



E ap (T c ) /T c > 45, 



(2) 
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where E sp and T c are the sphaleron energy and the crit- 
ical temperature, respectively. Since it was shown in Q 
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that Es p (T) oc v (T) 1 , the condition (J2]) can be trans- 
lated for the case of Standard Model to 

v c /T c > 1, (3) 

where v c is the field value at the critical temperature. 
However this condition ([3]) is not fulfilled in the case 
of Standard Model, because the thermal induced cubic 
term 2 is not large enough; also this leads to an unaccept- 
able upper bound on the Higgs mass [12] 

m h < 42 GeV. (4) 

The constraint gets even stronger when the two-loops 
effect and a proper treatment of the top quark are in- 
cluded [f3| ■ It is clear that this bound is in contradiction 
with the lower bound coming from LEP m/j > 114 GeV 
0. 

But this severe bound can be avoided by adding a com- 
plex scalar singlet that couples only to the Higgs doublet 
with an appropriate choice of the theory parameters in a 
way that the singlet does not develop a vev 

If a new scalar (or many scalars which can be singlets or 
in doublet w.r.t SU(2)l) acquiring a vacuum expectation 
value x are added to the Standard Model, the term v c 
in (J3J should perhaps be replaced by f2 c which equals 
{v 2 + x 2 } 2 ; or {v 2 + (x — 2£_) c } 5 when the false vacuum 
is [0,xq) instead of (0,0) [la] . Then Q becomes 

n c /T c > 1. (5) 

If the new particle (s) is a singlet (s), cubic terms can 
exist in the potential at tree-level, and therefore the phase 
transition gets stronger without the need of the thermally 
induced one 

In this work, we want to check for a model with a sin- 
glet whether the passage from ^ to © is true or not? 
We will consider the simplest extension of the Standard 
Model with a real singlet. This paper is organized as 
follows: In the second section, we introduce briefly this 
model, and find the sphaleron solution in the third sec- 
tion. In the fourth section, we discuss the strength of the 
first order electroweak phase transition (EWPT). And 
finally we give our conclusion. 



II. THE STANDARD MODEL WITH A 
SINGLET 'SM+S' 

Let us consider an extension of the Standard Model by 
a singlet real scalar S coupled only to the standard Higgs. 



1 This was also checked for the Minimal Supcrsymmctric Standard 
Model (MSSM) in [hJ, then (3) is valid also for the MSSM. 

2 This term is forbidden by symmetry at tree level, however it 
appears as T(m 2 )5 from the thermal bosonic contributions to 
the effective potential at one-loop. 



We concentrate here on the scalar sector (SM Higgs and 
the added singlet) and the SU(2)z gauge sector. 3 

The effective Lagrangian 

The Lagrangian is given by 

C = -\F* y F a ^ + {D^{D^) + \{d^S){d»S) 

-VefffaS), (6) 

where <j> is the Higgs doublet 

T = l/\/2 ( xi + % h + i X3 ) (7) 

where h is the scalar standard Higgs, x's are the three 
Goldstone bosons, and F^ v is the SU(2)l field strength 

F° v = d^Al - %Al + g<? hc A\Al. (8) 

Dfj, is the covariant derivative; when neglecting the Uy (1) 
gauge group, it is given by 

= a„ - ha° a K- ( 9 ) 

Finally V e f f (</>, S) is the effective potential, which is at 
tree-level given by 

V (<f>,S) = A|0| 4 - M 2 \<j 3 \ 2 + u:\^\ 2 S 2 + p\cp\ 2 S 

+^fS* 4 - fS 13 - ^fS 2 . (10) 

We can eliminate /^ 2 and by making (v,x) as the 
absolute minimum of the one-loop effective potential at 
zero temperature, where v = 246.22 GeV is the standard 
Higgs vev. 

Now, we write the explicit formula of the one-loop ef- 
fective potential. We will consider the contributions of 
the gauge bosons, the standard Higgs h, the singlet S, 
the Goldstone bosons Xi,2,3 and the top quark. The field- 
dependent masses at zero temperature are given by 

m' 2 - ln 2 h 2 m 2 - SL+Rlh 2 m 2 - alh 2 

to 2 = Xh 2 — /i 2 , + loS 2 + pS 
mis mh = |[(3A + w) h 2 + (3A S + u) S 2 
+ (p-2a)S- nl - fi 2 s 
T { ((3A -uj)h 2 ~ (3A S - u) S 2 + (p + 2a) S - p\ + p 2 s f 

1 

+A{2ujS + pfh 2 } 2 ] (ff) 

where yt is the Yukawa coupling for the top quark, and 
g 2 = g 2 + g 2 , however we neglected the t/y(l) gauge 
group and therefore g = and mz = mw ', and m 2 2 are 



Since we are interested here in the Sphaleron solution, we as- 
sume that the Uy ( I) contribution to the sphaleron energy to be 
negligible as in the case of Standard Model [lH . 
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the Higgs-singlet eigenmasses. Then the one-loop cor- 
rection to the effective potential at zero temperature is 
given by 



V^°(h,S) = niG(mUh,S)) (12) 

i=W,Z,t,h,S, X 

,2 r / „ \ 3- 



and G (x) 



Here Q is the renormalization scale, which we take to be 
the standard Higgs vev Q = v; and m are the particle 
degree of freedom; which are 

n w = 6, n z = 3, n h = 1, n x = 3, n s = 1, n t = -12. 

The temperature-dependent part at one loop is given 



V?*°{h,S,T) 



rp4 

2^2 



and 



Y ni JB,F(ml(h,S)/T 2 ) 

i=W,Z,t,h,S,x 

(14) 



J B (0) = J dx a; 2 log|l-exp -y/x 2 + 6 } 

- - 45 + 12° - 6° ~32 l0 S^' ( 15 ) 
poo 

J F {6) = J dx x 2 log |l + exp |- \/x 2 + 6>J j 

- 360 " 24^ - 32 l0 gZJ' ( 16 ) 

where a& = 167r 2 exp(3/2 — 2te), a/ = tt 2 exp(3/2 — 27^) 
and 7s = 0.5772156649 is the Euler constant. There is 
also another part of the effective potential which is the 
ring (or daisy) contribution [20) 

V ring (h,S,T) = nA{M 2 {h,S,T)) 3 * 



i=W,Z,h,S,x 



(17) 



where M 2 (h,S,T)'s are the thermal masses of the 
bosons, which are given by 



Mf (ft, S, T) = m\ (ft, S) + n !: 



(18) 



and Ilj is the thermal correction to the mass, its values 
for the bosons in our model are: 



Ll w — 


11 9 , 

~6~ 






u L - 




n| = o 




n x = 


(.9 2 /4 + 


A/2 + y 2 /4 -+ 


w/6) T 2 


Tlhh = 


(.9 2 /4 + 


3A/2 + y 2 /4 


+■ w/6) T 2 


n S 5 = 


(A s /4 + 


2w/3)T 2 , n 


^5-0 



where the script L (T) denotes the longitudinal (transver- 
sal) mode for W and Z. Then the full one-loop effective 
potential at finite temperature is the sum of (|10|). (| 1 2 [1 . 
(O and CE7]): 

K//'° OP (*> £ T) = ^0 (ft, 5) + ]T (m 2 (ft, S)) 

i=W,Z,t,h,S,x 
+ J2 n l J B ^ F (m 2 {KS)/T 2 ) 



2tt 2 



i=W,Z,t,h,S,x 



--!- g { (M 2 (ft, 5, T)) 1 - (m 2 (ft, 5)) 1 } . 

(20) 



i=W,Z,h,S, x 



The mass-squared values of the Goldstone bosons or 
the Higgs-singlet eigenstates can be negative. In the case 
where a mass value (or more) is negative, the cubic term 
in (|15[) becomes non analytic, that's not a problem since 
it is already replaced by the thermal mass in (| 1 T[> . where 
it will be compensated by the thermal correction. If the 
thermal correction is not enough to compensate the neg- 
ative value, this term should be omitted since it is imag- 
inary and does not belong to the effective potential. 

The space of parameters 

In our theory, we have quite a few parameters, 

A, As, to, p, a, p 2 h , p%, 

in addition to the singlet vev x. As mentioned above, p? h 
and p? s can be eliminated as 



/4 



= Xv + lux + px + — 



> 2 + £ £ + A s x 2 



dh 



h—v 
S=x 

dV?=°{h,S) 



as 



(21) 



a 22 ) 

S—x 



(19) 



after which our free parameters are A, As, ui, p, a and 
x. Since the theory is invariant under the discrete sym- 
metry (x,p,a)—>(-x,-p,-a), we will assume only positive 
values for the singlet vev x. We want also to keep the 
perturbativity of theory by imposing A, As, \u\ <C 1. We 
choose the parameters, A, As, u>, p, a and x, lying in 
the ranges: 

0.001 < A, As < 0.6 
-0.6 < to < 0.6 

100 < x/GeV < 350 
-350 < a/GeV < 350 
-350 < p/GeV < 350 (23) 

The stability of the theory implies that the potential 
goes to infinity when the field goes to the infinity in any 
direction, which implies to 2 < A x As- Moreover, we need 
that any minimum or extremum of the potential should 
be in the range of the electroweak theory; let us say that 
all the minima and extrema must be inside the circle h 2 + 
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S 2 = {600 GeV} 2 in the h— S plane; and therefore the 
potential is monotonically increasing outside this circle 
in any direction. 

In the Standard Model, the Higgs mass lower bound is 
given by m^ M > 114 GeV [14(. The mixing between the 
standard Higgs and the singlet changes the couplings of 
the standard Higgs to all the SM sector (gauge bosons 
and leptons), and therefore this bound is not viable. In 
our work, we will not derive the new lower bound for 
the Higgs mass, but we will restrict ourselves only with 
masses mi, 2 in the range 65 GeV to 450 GeV. 



III. SPHALERON IN THE 'SM+S' 

In order to find the sphaleron solution for this model, 
we follow the same steps as in the SU(2)l model. Ap- 
plying Euler-Lagrange conditions on the effective La- 
grangian, (fl"2")l or ([20)1 . we find the field equations 

d^F^ T - qe qab A b n F aaT + \q 2 h 2 A qT = 



d 2 h - \g 2 hA a ^ + ^Veff (0, S, T) = 
9 2 S+^V ff (<f>,S,T) = 0. 

We will work in the orthogonal gauge where 

A a = 0, XiAi = 0. 



(24) 



(25) 



We will not use the spherically symmetric ansatz for 
{(f>, Ai} in 7], but another equivalent one [2lj |. 



A? (a) = 2(l-/(r))-=^ 
gr z 

H{x) = ^(J). h = vL(r) 
S(x) = xR{r). (26) 



Here v and x are the Higgs and singlet vevs in the general 
case (zero or nonzero temperature) . Then one can rewrite 
the field equations (T2"4")l as 4 



C 2 ^/ = 2/(l-/)(l-2/)-I^C 2 i 2 (l-J) 



d <2 d 



2£(i-/r 



C 2 dVe ff (h,S,T) 



jLr 2 — R = £ av eff (h,s.T) 

dC d( g 2 xn 2 ds 



where £ = gflr; the parameter f2 can take any non- 
vanishing value of mass dimension one (for example v, 
x or \Jv 2 + x 2 ); and the energy functional is given by 



h—vL^S—xR 



h=vL,S=xR 

(27) 



4 There is a similar work done in [22| , however there is a difference 
in the definition in the theory parameters, and also there is an 
error on the r.h.s of the first equation in (19) in this paper, 
where the term u 2 /v 2 should be corrected as u 2 /V 2 according 
to his notation. In our notation it is the term v 2 /D. 2 in the first 
equation in 11271 . 



(i-/r 



v 2 » 9 / d \ v 2 



+ ^ C2 \dC R J ' g 2 ^ 
{V eff (vL, xR, T) - V eff (v, x, T)}} , (28) 

with the boundary conditions (See Appendix A) 



for C - 



L 



R 



c 2 

+ K 



for C 



00 



2. 



/- 
L - 
R- 



1 
1 
1. 

(29) 

Let us now compare the energy functional ()28|) to that 
of the minimal Standard Model (eq. (10) in [7j). The 
difference between these quantities is of course the con- 
tribution of the singlet, which contains the kinetic term, 
the mixing with the standard Higgs; and a contribution 
to the potential term. However if we compare ()28p with 
the same quantity in the MSSM case (Ea. (2.22) in [l(|), 
we find that in the MSSM both Higgs fields, hi and h 2 , 
have similar contributions to the sphaleron energy, and 
its general form remains invariant under hi <-*■ h%. How- 
ever this is not the case for ([28)1 if h S, because of a 
missing term like R 2 (1 — /) . 5 

For the MSSM sphaleron energy, its form is invariant 
under hi h 2l and it scales like {vf + v\ } 2 ; and for our 
model 'SM+S', a similar invariance is absent. Could it 
nevertheless be that E$ p oc {v 2 + (x — xq ) }5? We will 
check this in the next section. 

But when comparing (|28[) with the same quantity for 
the Next-to-Supersymmetric Standard Model (NMSSM); 
(eq. (2.20) in [23|; after eliminating explicit CP phases), 
we find no large difference expect for what comes from the 
fact that the NMSSM contains a doublet more than the 
'SM+S'; and we remark also similar equations of motion 
and also similar boundary conditions. 

The analytic solution of the system ([2"7| is not possible, 
this should be done numerically. To solve this system 
numerically, we need to transform it into a system of 6 
first-order differential equations, and therefore we have a 
first order two-point boundary problem, then we use the 
so-called relaxation method to solve it. This method is 
well explained in section 17.3 of (24|. 

As an example, we solve the system ([2T)l for four chosen 
sets of parameters (A, B, C and D); and then we can 



To be more precise, the absence of a mixing between the singlet 
and the gauge field is not the only reason to spoil this invariance, 
but this invariance is absent also in the tree-level potential. 



5 



compute the sphaleron energy (|28|) at any temperature 
T < T c . All the results for the sets A, B, C and D are 
summarized in table |H 
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A 


B 


C 


D 


A 


0.4000 


0.4000 


0.5000 


0.4150 


As 


0.4003 


0.4200 


0.4100 


0.5500 




0.3818 


0.2818 


0.3818 


0.3000 


x/GeV 


200 


250 


350 


350 


a/GeV 


-38.89 


38.89 


38.89 


194.44 


p/GeV 


-272.22 


-194.44 


-272.22 


-300 


mi/GeV 


178.00 


204.00 


244.74 


203.05 


m 2 /GeV 


311.92 


269.80 


333.96 


318.80 


T c /GeV 


141.55 


241.34 


389.94 


270.08 


E sp (0)/GeV 


9618.6 


9721.3 


9883.3 


9726.6 


Vc/T c 


1.680 


0.838 


0.495 


0.386 


n c /T c 


3.138 


1.232 


1.436 


0.703 


E Sp (T c )/T c 


64.851 


32.980 


20.459 


13.577 



Table I: Representative parameter values and the correspond- 
ing values of the scalar masses, critical temperature and dif- 
ferent ratios needed for the criterion of a strong first order 
phase transition. 



From table [H the set (A) satisfies both conditions (O 
and |(3J), (D) does not satisfy either of them, and both 
(B) and (C) satisfy © but not ©. 

The profiles of the functions /, L and R are shown in 

Fig. rjj 




T (GeV) 



r (GeV) 



Figure 1: A, B, C and D represent the profiles of the functions 
f, L and R for the sets of parameters A, B, C and D in table 
[J respectively. The continuous lines represent the profiles at 
zero temperature and the dashed ones represent the profiles of 
the functions at finite temperature. 

From all cases in Fig. [TJ we remark that the singlet 
profile is not much different than the unity, due to Neu- 
mann type boundary at r = 0. This comes from the fact 



that the singlet couples to the Higgs doublet and not the 
gauge fields. Then we claim that the singlet contribution 
to the sphaleron energy (|28p . should be small compared 
to doublet and gauge field contributions. 



IV. THE PHASE TRANSITION IN THE 'SM+S' 

In Ref [3, the authors have studied the EWPT 
strength using the same tree-level potential as (|10l) with 
some differences in the parameter definitions. They eas- 
ily got a strong first order phase transition even for Higgs 
masses much larger than (f4|). And of course they used the 
criterion ([5]) instead of ([3]), where the quantity v c is re- 
placed by f2 c = {vl + (x c — x 0c ) }5 . Since Vl c /T c > v c /T c 
is always fulfilled, the phase transition gets stronger for a 
larger parameter space compared with the minimal Stan- 
dard Model case. 

Let us take a random choice of about 3000 parameters 
in the ranges (|23p . and make a comparison between the 
two different criteria of the strong first order phase tran- 
sition ([5]) and (J2j) . We show the plots of the quantities 
il c /T c and E$ p (T c ) /T c as functions of the lightest Higgs 
mass mi in Fig. O 

Comparing the number of points above and below the 
dash-dotted line in both cases (a) and (b) in Fig. [21 we 
remark that the first order phase transition is stronger 
than that of the Standard Model with both criteria. How- 
ever according to the large number of points below the 
dash-dotted in (a), there are a lot of points which satisfy 
([5|) but they do not really give a strong first order phase 
transition according to @ . 

When comparing the points in Fig. [2}-a with the curve 
which represents the Standard Model case, we remark 
that the addition of a singlet increases, in general, the 
quantity Es p (T c )/T c which is relevant to the phase tran- 
sition strength; that there are even a large number of 
points above the line Es p {T c )/T c — 45. 

The passage from the criterion j2]), which is model- 
independent, to ([31), was based on two assumptions [TT[ | : 

(I) The sphaleron energy Es p {T) scales like the vev 
v(T). 6 

(II) The sphaleron energy at T = 0, is taken to be 1.87 
in units of Anv/g. 

If the assumption (I) is satisfied in our model 'SM+S', 
i.e. E Sp (T) oc fi(T); and E Sp (0) ~ 1.87 x 47rO (0) /g, 
then ([5]) is the condition of a strong first order phase 
transition, but this not the case as mentioned above. 

In general, the value of the sphaleron energy at zero 
temperature is significantly different from 1.87 in units of 
(4ir£l (0) /g), thus if the assumption (I) is fulfilled, then 
the criterion (HI) is still viable but should be relaxed as 



6 As mentioned above, this was verified for the Standard Model 
and the Minimal Supersymmetric Standard Model 10]. 




50 100 ISO 200 250 300 




(b) 



50 100 150 200 250 300 



Figure 2: For the points above the dash-dotted lines in (a) 
and (b), the electroweak phase transition is strongly first order 
according to (TJ|) and (5$), respectively. In (a), the continu- 
ous curve represents Esp(T c )/T as a function of the Higgs 
mass for the case of the Standard Model. 

^c/T c > 1 + 5, where 6 describes the deviation from the 
assumption (II). 

In order to probe the assumption (I) for our case, i.e. 

E Sp (T) oc Q (T) , (30) 

we take the sets (A), (B), (C) and (D) used in table U 
in the previous section, and plot the ratios v (T) jv (0), 
(T) /Q (0) and E Sp (T) /E Sp (0); as functions of tem- 
perature, which lies between the critical temperature and 
another value. The results are shown in Fig. [3l 

Let us here comment on Fig. [3] For the case of (A), 
the ratio E Sp (T)/ E Sp (0) is close to both v(T)/v(0) and 
f2(T)/O(0), which is almost 1 at T c . For the case of 
(B), the ratio E Sp {T)/ E Sp (0) is very close to v{T)/v(0); 
however it is different a bit from 1 at T c . At the temper- 
ature T ~ 204.5 GeV, there exists a secondary first order 
phase transition, it happens on the axis h = 0, where the 
false vacuum (0, Xo) is changed suddenly. In the case (C), 
the ratio E Sp (T)/E Sp (0) is closer to v(T)/v(0) than to 
n(T)/ft(0); it is also different from 1 at T c . In the last 
case (D), there is also a secondary first order phase transi- 
tion around T ~ 256 GeV; where the true vacuum (v, x) 
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320 330 340 350 360 370 3S0 390 200 210 220 230 240 250 260 270 



T(GeV) T (GeV) 

Figure 3: The solid line denotes the ratio v(T)/v, the 
dashed one denotes Q (T) /fi; and the dot-dashed one denotes 
Esp{T)/ Es P (0) ■ All the plots end at the critical temperature. 



changes discontinuously. We cannot call this an EWPT 
because the scalar h has already developed its vev. The 
ratio E Sp (T)/E Sp (0) is still scaling like v(T)/v(0), but 
significantly different from 1 at T c . 

It is clear that Es p (T) does not scale like fi(T), but 
roughly speaking it scales like v(T); with a little devia- 
tion in some cases. 

We claimed previously that the contribution of the sin- 
glet S to the sphaleron energy is small, and therefore 
this may be the reason why Es p {T) does not behave like 
f2(T); and also does not behave exactly like v(T). In 
order to estimate the effect of the singlet field S on the 
sphaleron energy (|28[) , we compute the sphaleron energy 
(|28p with replacing the singlet S by its vev x, which we 
denote £s P (T). Then the fifth term in (f2"5| disappears 
and the third equation in (f2"4l and [2T|) disappears also; 
the problem is reduced to a SU(2)l Higgs-gauge system 
Q; but with a modified potential V e ff(h)=V e ff(h,x,T). 
We find that the singlet S gives a negative contribution 
to the sphaleron energy which is larger at higher temper- 
atures. But the contribution size is, generally, negligi- 
ble. The maximum contribution in the case (A) is less or 
equals 2.1 %, less than 1.1 % for in case (B), and almost 
zero in the case (C): it is less than 0.08 %; and this is 
expected because the function R in Fig. [T]-C is very close 
to 1. In the case (D), there are two different phases: at 
the first one before the secondary phase transition i.e, 
T c > T > 256 GeV; the singlet contribution is significant 
(between 8~9 %), this may be due to the smallness of 
the Higgs doublet vev at this range. While in the sec- 
ond phase T < 256 GeV; the singlet contribution, as in 
the other cases, is less than 2 %. Then in the absence 
of secondary first order phase transitions, we can neglect 
the singlet contribution, but in its presence the singlet 
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contribution can be sizeable but not as large as that of 
the Higgs doublet or gauge fields. 

To justify this picture, we take again 3000 random sets 
of parameters and plot Es p (T c ) /T c as a function of Q. c /T c 
in Fig. |3 
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Figure 4: Es P (T c )/T c vs Q. c /T c for 3000 randomly chosen sets 
of parameters. 

Since there exist too many points in the region 
(E Sp (T c )/T c < 45 n fl(T c )/T c > 1), the criterion © is 
not the definition of a strong first order EWPT. However 
it is satisfied for all points that give really a strong first 
order EWPT except for 10 points due to the existence 
of secondary first order phase transitions. Then we are 
now sure that ((3J) does not describe a strong first order 
EWPT. 

In the sphaleron transitions, the singlet S has no rela- 
tion to lepton or baryon number breaking phenomena. It 
does not couple to fermions or gauge bosons; it is just a 
compensating field in the field equations; l|24p and (|27p: 
and its effect on the sphaleron transition is negligible as 
shown above. Then we claim that only the Higgs doublet 
vev is relevant for the phase transition strength. 

We take 3000 random sets of parameters used previ- 
ously, and plot Es p (T c )/T c as a function olv c /T c in Fig. 
® 

It is clear that Es p {T c )/T c scales almost exactly 7 
like v c /T c except for some points, and © can describe 
the strong first order EWPT criterion for most of the 
points. Then when studying the EWPT in models with 
a gauge singlet, one should treat the problem as the 
SM case (in case of one doublet) with replacing the sin- 
glet by its vev; and look for the Higgs vev in the path 
dV e ff(h,S)/dS\ s=x — 0; whether it is larger than the 
critical temperature i.e, v c /T c > 1? 



Except some points due to the existence of secondary first order 
phase transitions; or due to the significant singlet contribution 
to the sphaleron energy; especially for smaller Higgs vev values. 
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Figure 5: Es p (T c )/T c vs v c /T c for 3000 randomly chosen sets 
of parameters. 

With this modified potential V e //(/i) = 
V e ff(h,S)\ s _ x ; the EWPT can be obtained easily 
as done by the authors in [251 ] . 

V. CONCLUSION 

In this paper, the electroweak phase transition for 
the Standard Model with a singlet is studied using the 
known criteria in the literature in addition to the model- 
independent criterion found in H . The authors [l(| E3] 
found that the EWPT gets stronger even for Higgs 
masses above the bound Jl]). They modified the sim- 
ple criterion ([3|) into (|5|), where they replaced the Higgs 
vev by the distance between the two degenerate minima 
in the h — S plan. 

In our work, we checked whether this criterion is viable 
for this kind of models or not. We took the Standard 
Model with a real singlet, then we studied the EWPT 
using the sphaleron configuration at the critical temper- 
ature, then we checked whether all the steps of the pas- 
sage from the model independent criterion (|2"|) to (|3"j) in 
the Standard Model case, are respected for our model 
(i.e. the passage from © to (JSJ) or not? 

We found that the EWPT gets stronger even for Higgs 
masses larger than 100 GeV; and this model does not 
suffer from the severe Higgs mass bound (j4j. However, 
we remarked that a sizeable number of the parameters 
sets satisfy the modified criterion but do not really give a 
strong first order EWPT, this allowed us to conclude that 
the 'modified criterion' is not the criterion that describes 
a strong first order EWPT. 

In order to understand why this modified criterion is 
not viable in this case, we returned back to the SM to 
see how the passage from the model-independent crite- 
rion to the simpler one proceeds? We found that the two 
assumptions needed for the passage to the simpler crite- 
rion are not fulfilled, in general, in our model 'SM+S': 

• The sphaleron energy at zero temperature is different 
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from the value 1.87 in units of (4nfl/g). 

• The sphaleron energy at finite temperature does not 
scale like O(T). 

We guess that the reason for this is that the sin- 
glet does not couple to the gauge field, then the miss- 
ing of some contributions to the sphaleron energy like 
R 2 (1 - if m (EHD, can spoil the scaling law, E Sp (T) cc 

(T). This can be inspired if we compare this situation 
with the case of the MSSM, where this scaling law does 
work; and the general form of the sphaleron energy is 
invariant under hi <-> h^. The fact that the singlet does 
couple only to the Higgs doublet leads to a singlet pro- 
file in the sphaleron configuration in a Neumann type at 
r = 0, which makes the singlet contribution too small. 
Another important remark is that the possibility of sec- 
ondary first order phase transitions can, sometimes, spoil 
this scaling law. 

As a conclusion, we can say that the condition Q c /T c > 
1 is not valid as a strongly first order phase transition 
criterion. But the usual condition v c /T c > 1, is still 
the viable one, which can describe the strong first or- 
der phase transition for the majority of the physically 
allowed parameters as stated in Fig. [5] Moreover, this 
can be satisfied even for Higgs masses in excess of 100 
GeV unlike in the Standard Model. 

Then in such a model where the singlets couple only to 
the Higgs doublets, it is convenient to study the EWPT 
within an effective model that contains only doublets, 
where the singlets are replaced by their vev's. We expect 
similar conclusion for models like the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM), where in 
this model the singlet couples only to the two Higgs dou- 
blets; and its profile is a Neumann type in the sphaleron 
configuration. Then the criterion for a strong first order 

EWPT is {v\ + v\Y l T > 1 at tne critical temperature, 

instead of {vj +v% + (x- x ) 2 } 3 /T > 1. 
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Appendix A: THE BOUNDARY CONDITIONS 



To find the boundary conditions of (|2T[) . one should 
take into account that the energy functional (|2"5)) should 



be finite. It is clear that in order for the contributions 
of the second and fourth terms in (J2HJ) to be finite, / 
must go to unity at the limit ( — > oo. According to the 
sphaleron definition, scalars go to their vacuum at the 
infinity, i.e. L, R — > 1 when ( — > oo, which makes the 
last term contribution to l[28p finite. Thus one can write 
all the functions as 1 — a exp {— c?iC}i an d find the values 
of Ci and di by inserting this behavior into the differential 
equations (|27|) . 

In the limit £ ~ 0, let us assume that the functions /, 
L and R have the profiles 



/(C) ~ c n/ 

L(0 ~ c 1+ C L 

R(0 ~ c 2 + C R , 



(Al) 



where rif, and tir are some positive constants. In this 
limit, ([27|) can be approximated as 



d( 2 
d( 2 



2/ 



L 



R 



2 / 



1 u T 2 



'm L ' 

2 d 



R 



dVeff(h,S } T) 

dS 



(A2) 



S=xR 



From the second equation in (|A2j) . one can easily con- 
clude that L ~ £ or ~ C -2 , however the second choice 
makes the energy functional integral (|28|) divergent, thus 
L ~ £ or { c i — 0? n L — 1}. Using this result, one can 
conclude from first equation in (|A2|) that / ~ ( 2 . How- 
ever the situation is different for the last equation (|A2[) . 
then one can make 



1 dV cff (h,S,T) 
g 2 xtl 2 dS 



h—vL 
S=xR 



a( 2 + {A + b( 2 }R(() 
+BR 2 (C) + CR 3 (C) (A3) 



then inserting (|A1[) in (|A3|) . one finds that the only pos- 
sibilities are tir = — 1 and tir — 2, where the first choice 
is excluded in order that the energy functional integral 
(|28p to be convergent, thus R ~ a + bC, 2 . Therefore at 
C = 0, R satisfies the boundary condition of Neumann 
type, while / and L satisfy those of Dirichlet type. The 
boundary conditions are summarized in (|29[) . 
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